Global minimizers for axisymmetric multiphase membranes. 

Rustum Choksi* Marco Morandotti^ Marco Veneroni* 

March 3, 2013 



Abstract 

We consider a Canham-Hclfrich-type variational problem denned over closed sur- 
faces enclosing a fixed volume and having fixed surface area. The problem models 
the shape of multiphase biomcmbrancs. It consists of minimizing the sum of the 
Canham-Hclfrich energy, in which the bending rigidities and spontaneous curvatures 
are now phase-dependent, and a line tension penalization for the phase interfaces. 
By restricting attention to axisymmetric surfaces and phase distributions, we extend 
our previous results for a single phase [7] and prove existence of a global minimizer. 

Keywords: Helfrich functional, biomembranes, global minimizers, axisymmetric sur- 
faces, multicomponent vesicle. 

AMS subject classification: 49Q10, 49J45 (58E99, 53C80, 92C10). 

1 Introduction and main result 

Lipid bilayers are the most elementary and indispensable structural component of bio- 
logical membranes, which form the boundary of all cells in living systems. In biological 
membranes, the bilayer consists of many different lipids and other amphiphiles, which 
can separate into coexisting liquid phases, or domains, with distinct compositions. Two 
types of phases are typically observed: a liquid phase with short-range order and a liquid 
disordered phase (see, e.g., [18], [4]), which we label phase A and B in the sequel. Their 
configurations, however, are fundamentally distinct from other interfaces, since they are 
not determined by a surface tension but rather by a bending elasticity, as introduced 
independently by [6], [15], and [13]. 

As in, e.g., [18] and [4], the Canham-Helfrich- Evans energy functional of a surface £ 
and a phase ip is given by 

^(S, <p) = (^-(H - H (<p)) 2 + nab) K^j dS + an\T), (1.1) 
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where dS denotes the integration with respect to the ordinary two-dimensional area 
measure, H is the sum of the principal curvatures of S, i.e., twice the mean curvature 
and K is the Gaussian curvature. In a state where the phases are completely separated by 
a sharp interface, we let 99 : S — > {0, 1} denote the characteristic function of the phase 
A, kh ((f), k-g (v 9 ) are the phase-dependent bending rigidities, and Ho(tp) is the phase- 
dependent spontaneous curvature. In the last term of (1.1), > is the (constant) 
line tension coefficient, T is the interface between the phases A and B, and Ti 1 is the 
one- dimensional Hausdorff measure. 
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Figure 1: Two examples of axisymmetric multiphase surfaces. 

We prove the existence of a minimizer of (1.1) in the set of families of axisymmetric 
surface-phase couples (£,</?) such that the total area of the surface S, the area ratio 
between the phases, and the volume enclosed by the surface E are fixed. 

In the last forty years, homogeneous membranes have been extensively studied from 
the experimental, theoretical, and numerical points of view (see, e.g., [21]). Investigation 
of inhomogeneous systems started more recently ([17], [18]) and is nowadays at the 
center of an increasing focus. Multicomponent vesicles have recently been considered 
in numerical studies aiming at understanding the dynamics of the phase separation, the 
stability of nanodomains, and the complex morphology of the membranes; see the review 
paper [12]. We refer in particular to the paper [22] concerning the dynamics of one- 
dimensional curves in viscous fluids, while two-dimensional surfaces where investigated 
with the phase field method in [24] and [11], by means of surface finite elements in [10], 
and by adaptive finite elements in [19]. Our result provides a theoretical basis, at least in 
the axisymmetric case, for the existence of the shapes approximated by these numerical 
studies. More advanced models, with respect to (1.1), which couple the chemical and 
mechanical properties of lipid bilayers, are subject of current research, see, e.g., [26], [8]. 

We note that in the case where the spontaneous curvature vanishes, the Canham- 
Helfrich functional becomes the famous Willmore functional (see, e.g., [25, Chapter 7]). 
While there has recently been tremendous amount of research associated with minimiz- 
ing the Willmore functional, this research does not directly carry over to the doubly- 
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constrained Canham-Helfrich functional (cf. [7]). 



1.1 Axisymmetric multicomponent vesicles 

We detail now the assumptions on the variables, the parameters, and the constraints 
that compose our problem and we state the main result. 



i) Surfaces. Axisymmetric surfaces, i.e., surfaces of revolution, can be obtained by 
rotating a curve about a line. Let IR+ := {(x,y,z) £ M 3 : x > 0, y = 0} be a half-plane 
in M 3 ; we first consider a smooth curve 7 : [0,1] (71(f), 0,72(t)). By rotating 

7 around the axis x = y = 0, we obtain the surface E parametrized by: 

r(t, 9) = [ 7l (t) cos(0), 71 (*) sin(0), 72 (t)] , (*, 0) G [0, 1] x [0, 2vr]. (1.2) 

If a surface E admits the parametrization (1.2), then we say that E is generated by 
7. A standard computation (see Section 2.2 below) shows that if 7 generates a smooth 
surface E without boundary, then the 2-dimensional surface area |E|, the enclosed volume 
Vol (E), and the principal curvatures fei, ki of the generated surface are given by 



2tt / 7i|7|dt, Vo1(E) = tt/ jfadt, 

Jo Jo 



(1.3) 



, 7271 ~ 7172) , 72 

h = TTTo , k 2 = ——. (1.4) 

m 3 71I7I 

Denoting by C 1 the one-dimensional Lebesgue measure, let /it 7 := 27T7i|7|£ 1 l[ 0j i] be the 
area measure induced by the curve 7. It is not difficult to see, at least in the case of 
only one phase (see [7]), that a bound on ^(E) provides an a priori estimate on 7 in 
the space L 2 ((0, 1); |U 7 ), which translates into a bound for 7 in the space W 2,2 on any 
stretch of curve such that 71 > e > 0. Precisely, by Proposition 3.2, it is not restrictive 
to assume that 

7GC 1 ((0,1);M2) ) 7i>q vte (0,1), 7ei 2 ((0,l);4 (1-5) 

Moreover, in order for ^(E) to be finite, 7 is forced to meet the z— axis orthogonally, 
so if 71(0) = 71(1) = 0, then r is the CA immersion of a closed surface. If r is not an 
embedding, as, for example, in the case of the surface generated by the limit curve in 
Figure 2, we define the generalized two-dimensional surface area, enclosed volume and 
principal curvatures of the generated surface by the quantities in (1.3)-(1.4). 



ii) Phases. If <p is the characteristic function of phase A, then in order for J^(E, (p) 
to be bounded, and in particular in order for aT-L l {T) to be bounded, we need to impose 
some kind of regularity on the class of admissible phases <p : E — > {0, 1}, for example 
by restricting to characteristic functions of finite perimeter sets on the surface E. Under 
the simplifying assumption of axisymmetry for the phases, as well as for the surfaces, 
a useful approach is then to follow the parametrization = ^07 : [0, 1] — > {0, 1}. 
Let J(<p) C (0, 1) be the set of points where if has a jump discontinuity. Owing to 
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axisymmetry, the (measure-theoretical) interface T between the two phases on the surface 
S is a union of circles 

r = {r(t,6) : t G J{ip), 9 G [0,2vr]} 

and 

U\Y) = 2tt! 7i(t)d||^||=27r ^ 71W, (1-6) 

where ||-Dv?|| = 'H lj( v ,) is the counting measure restricted to the jump set of 99 (see 
Section 2.3 below). Since 71 G C°([0, 1]) and 71 (t) > for t G (0, 1), then according to 
(1.6), ^(r) < +00 if and only if 71 G L 1 (||D( / ?||). In the sequel, it will be convenient to 
deal with the weaker request that 

<p G BV\ oc ((0, 1); {0, 1}), (1.7) 

extending & to +00 if the quantity in (1.6) is not bounded. The area measure of phase 
A can then be expressed as 

f (pdS = 2ir [ <p(t)7i(t)|7|dt. (1.8) 
Jo 

It would also be possible to choose, as ambient space for the phases, the space of special 
functions of (locally) bounded variation SBV\ oc . Since we are dealing with two-valued 
functions, and the Cantorian part of the measure will not appear in any case, we prefer 
to keep the setting as simple as possible and use BV\ oc functions. 



iii) Parameters. Let the bending rigidities for phase i be given by n l H , k g , and the 
spontaneous curvature be Hq , for i = A, B. In order for the the functional & to be 
coercive, we require (see Lemma 3.1 below) that 

k 1 h >0 and -p€(-2,0) for i = A,B. (1.9) 

K H 

We note that the physical range in which the parameters kh and kq are typically found 
is contained in the one we impose in (1.9), see e.g. [23] and [3] (note that the latter cites 
the former, but inverting numerator and denominator, by mistake). 

We define the coefficient functions on the interval [0, 1] in such a way that kh is the 
linear interpolation of kh{0) = 1 k h(1) = k-h > an d t ne same holds for kg and Hq , 
i.e., 

K H (<p):=<pK%+(l-<p)K%, kg(<p)--=<Pi4 + 0--<p)"G, %) :=^ + (l-# B 

According to the parametrization (1.2), the contributions in (1.1) depending on the 
mean curvature, the Gaussian curvature, and the line tension of the interface between 
the phases, can then be written as 

^ J KH (<f)(H - H ( V )) 2 dS = tt^ 1 k h (<p) (fci + k 2 - H ((p)) 2 7l | 7 | dt, (1.10) 
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f KG(<p)KdS = 2tt [ K G (<p)k 1 k 2 'ri\'y\dt, (1.11) 
Jt, Jo 

aH 1 {T) = 2aTT [ -fi(t)d\\D^\\. (1.12) 
•/ o 

Definition 1.1. We say that a couple of surface and phase (£,¥>) is admissible if the 
surface X is generated by a curve 7 satisfying (1.5) and the phase ip satisfies (1.7). 

Finally, for fixed area and volume, a configuration made of several connected com- 
ponents may have a lower energy than a one-component configuration. This could also 
be favored by a relatively high value of a, since separation of phases in different com- 
ponents would have no interface between phases and thus aTHMV) = 0. From a dy- 
namical point of view, in certain conditions, shape transformations involving topological 
changes, like budding and fission (see, e.g., [21, Section 3]), could be expected. We 
take this possibility into account by studying families of admissible surface-phase cou- 
ples S = {(Si, fi), . . . , (S m , (fm)}, and defining the total energy of such a system as the 
sum of the Helfrich energies of the single components: T(S) := YliLi ^"(Sj, Vi)- 

iv) Main result. We prove the following. 

Theorem 1.1. Let n l H , k g G R, i = A, B be given such that (1.9) is satisfied. Let a > 0, 
H'q £ R, i = A, B be given. Let A, IFa, V > be given such that 

J3/2 

V<—=, 0<n A <A (1.13) 
by 7r 

Let S(A, Ha, V) denote the set of finite families S = {(Si, 931), . . . , (Y< m ,{p m )} of admis- 
sible couples of surfaces and phases such that the generalized area, volume, and phase 
area constraints 

J2\^i\ = A, ^Tvoi^)^ X> / <p i (t)r n (t)\j\dt = n A 
1=1 1=1 i=i ^° 

m 

are satisfied (see (1.3), (1.8)). Let F(S) := ^"(S^, <pj). Then the problem 



i=i 



min{^(,S) : S e S(A,U A ,V)} 

has a solution. 



1.2 Structure of the proof and plan of the paper 

The proof of Theorem 1.1 follows the direct method of the calculus of variations: given 
a minimizing sequence of (systems of) surfaces and phases S n satisfying the area and 
volume constraints, by compactness we obtain a subsequence converging to a system S, 
and by lower semicontinuity of the functional T we prove that S is a global minimizer. 

Regarding compactness, it is fundamental that the phase-dependent parameters kh , 
kg , Hq are chosen in such a way that the functional & is an upper bound for the L 2 -norm 
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of the second fundamental form of the surface S (Lemma 3.1). In this way we can exploit 
the compactness result obtained in [7] in the case of homogeneous membranes. The 
main idea in modeling the phases (see also Section 1.3) is to follow the parametrizations 
7 instead of their images, so that the sequence of parametrized phases is defined on a 
fixed interval rather than a sequence of surfaces. The drawback with this approach is 
that the phases are only of locally bounded variation; where the curves touch the axis of 
revolution, the area measure vanishes. In other words, where the horizontal component 
71 of the curve 7 becomes infinitesimal, the line tension part of the functional (1.12) may 
allow for an infinite number of discontinuities in the phases ip. However, this will not 
constitute a problem in the proof of lower semicontinuity, since the combined phases and 
curves are well-behaved (Lemma 3.5). 

The proof of the lower semicontinuity for the curvature terms (1.10) and (1.11) re- 
quires a special care, since we have to pass to the limit simultaneously in the surfaces and 
in the phases defined on the surfaces. A useful tool can be found in the function-measure 
pairs introduced in [16]. 

In Section 2 we describe the notation, we derive the geometrical quantities involved 
in Helfrich's functional, and we recall the basic definitions and the main results regarding 
functions of bounded variation and measure- function couples. In Section 3 we study the 
compactness and lower semicontinuity of a bounded admissible sequence, and end with 
the proof of Theorem 1.1. 

1.3 Discussion 

Modeling the phases. The choice of modeling the phases by following the parametriza- 
tion ip = (p o 7 : [0, 1] — > {0, 1}, instead that by functions with support on the surface 
{0, 1}, yields first of all a much simpler setting, since the domain becomes a 
fixed real interval, instead of surface. It is also a way to solve a more intrinsic question 
related to the possible ill-posedness of the problem. Consider, as in the example in Figure 
2, a sequence of curves 7 n and phases ip n such that, in the limit, two stretches of curve 
carrying different phases overlap. How should the limit phase be defined in this case? 




—a a x 

Figure 2: The problem of defining the phase on a curve when overlapping can arise as 
limit of well-defined configurations. 
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Any {0, l}-valued phase <p denned on the support of 7 would not be a limit (not even in 
the sense of distributions) of the sequence tp n . Instead, defining the parametrized phases 
on the interval [0,1], the situation in Figure 2 could be described, e.g., by choosing a 
constant sequence ip n equal to the characteristic function of the interval (0, 1/2), which 
clearly converges to the same characteristic function. 

Another way to approach the problem of overlapping curves could be, for example, 
to consider the varifolds [16] associated with the curves and to describe the segment 
[—a, a] x {0} by a stretch with double density. In the case of a homogeneous curve in two 
dimensions, this approach was followed, for example, in [5], but in the case of multiple 
phases it is not evident which approach to follow. A possibility could be by multivalued 
functions on varifolds, but we deem more natural to follow the parametrization of the 
curves. In higher dimensions this could be generalized by following the immersion of a 
Riemann surface as opposed to studying its image. 

Generalizations - surfaces. Experimental evidence of higher-genus membranes is 
known for homogeneous membranes (see, e.g., [21, pages 19, 20]), but we are not aware 
of any experiment or simulation showing multiphase biological or artificial membrane 
of genus higher than zero. Nonetheless, from a mathematical viewpoint, it is straight- 
forward to allow for genus- 1 axisymmetric membranes (see [7]) by considering closed 
generating curves with strictly positive first component. Since the W 2 ' 2 -norm of such 
curves is bounded from above by (1.1), they have the same compactness properties of 
the curves that we consider, and Theorem 1.1 could be directly extended to include also 
genus-1 generators. What is not clear to us, is whether the genus of the minimizer could 
be prescribed, since a priori a sequence of genus-1 minimizer may degenerate to a genus-0 
surface in the limit. 

Generalizations - phases. Even though, according to experiments and simulations, 
we expect the phases of an axisymmetric surface to be axisymmetric as well, one could 
pose the problem of studying general phases on surfaces of revolution. We believe that 
this step could be performed following the same steps as in the symmetric the 
theorems used for compactness and semicontinuity for functions of bounded variation do 
not depend on the dimension. Actually, the functional (1.1) for a two-dimensional phase 
would provide a better bound, namely in BV instead of BV\ oc . 

Generalizations to non axisymmetric surfaces and phases would require a completely 
different approach: our method relies on Ascoli-Arzela compactness for equicontinuous 
curves, which in one dimension can be applied owing to the compact immersion W 1 ' 00 
W 2 ' 2 , which fails in higher dimensions. 

Necessity of constraints. In order to model realistic configurations of multiphase 
membranes, we considered classes of surfaces with fixed total area, phases area, and 
enclosed volume. We notice, though, that the only constraint which is necessary in 
order to obtain compactness of a minimizing sequence is that on the total area, which 
is needed in order to bound from above the full second fundamental form of the surface 
(Lemma 3.1) and to ensure that the limit of a minimizing sequence does not vanish. The 
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area of each phase is then bounded by the total area, and the enclosed volume results 
bounded by isoperimetric inequality. Theorem 1.1 could then be restated by saying that 
if the total area is bounded from above and away from zero, then there is a minimizer 
of (1.1), and, in particular, there is one for any admissible choice of volume and phase 
area constraints. 



2 Preliminaries 
2.1 Notation 

We refer to [1] for the definitions of the following objects. Let C d and H d be the d- 
dimensional Lebesgue and Hausdorff measures, let C d \-A an d 1-L d ^A be the restrictions to 
the set A of the (i-dimensional Lebesgue and Hausdorff measures, respectively. For an 
open set Bel, let C C (E) be the space of continuous functions with compact support 
in E, let Cq(E) be the closure of C C (E) with respect to uniform convergence, and let 
RM(E) be the space of Radon measures on E, which can be identified with the dual of 
C C (E). If /j, € RM{E) and fi[E) < oo, we say that fi is a finite Radon measure. For 
a curve 7 : [0,1] — > M?, t t— >■ (71 (t), 72(^)5 we use the shorthand notation {71 > 0} to 
denote the set {t £ [0, 1] : 71 (t) > 0}, we denote the first and second derivatives by 7,7, 
and we define the measure ^ 7 G RM(R) by /i 7 := 27T7i|7|£ 1 l[ 01 ]. 



2.2 Parametrization 

The derivation of the principal curvatures for a surface of revolution can be found, e.g., 
in [9, Section 3-3, Example 4] or in [7]. Recalling the parametrization introduced in (1.2), 
we define the tangents to the surface along the coordinate lines as 



d_ 

dt 




r(t,9)= [7i(t)cos0, 7i(*)sin0, 72 (t)], 



r e :=— r(t,0)= [-71 (t) sin 0, 71 (t) cos 5, 0] , 

and notice that r t ■ re = 0, i.e., they are always orthogonal to each other. The first 
fundamental form is given by 



9(t,0) = 



' E F ' 




rt -n r t - r e 




' \m\ 2 





F G 




r g -r t re-re _ 







71 (t) 2 . 



sfg:= A /det(^)=7i(t)| 7 (i)|. 



The normal vector to the surface can be oriented either inwards or outwards, depending 
on the direction of 7. For curves parametrized in counterclockwise direction, it is inwards: 



n(t,0) 



rt x r e 
y/9 



■ [ - 7i(*)72(i) cos0, -7i(*)72(i) sinO, 7i(t) 7 i(t)] 



7i(*)l7(*)l 

pr^yj [-72 (t) cos 0,-72 (t) sin 6, 7l (t)]. 
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Now we start to compute the second fundamental form, using a constant-speed parametriza- 
tion. 



d 1 

n t : = ^ t n ( t ' 6 ') = j^yj [-72(i)cos6>, -7 2 (t)sin6», 71 (t)] 

d 1 
n e := ^ n (*'^) = TTT^T [72W sin6», -7 2 (t) cos0, 0] . 



W^'-' | 7 (t)l 
The second fundamental form is given by 



L 


M ' 




' n t 


n 


n t ■ r e 


1 


7271 - 7i72 





_ M 


N _ 




n e 


n 


n e -r e _ 


|7(*)I 





7i72 



II(t,6) 

The signed curvature of 7 is defined as 

k := 



727i - 7i72 



l7l J 



The Gaussian curvature is given by 

LN - M 2 (7172) (7271 - 7i72 ) 72^ 



K = hk 2 



EG - F 2 



(7i) 2 l7l 4 



7il7l 



The mean curvature is (for sake of notation, we define as mean curvature the double of 
what is often defined as mean curvature) 

, _ LG - 2MF + NE _ 7^(7271 ~ 7172) + 7172I7P _ , , 72 
^" fel + fc2 " EG — F 2 ' ( 7l )2| 7 |3 + 

The principal curvatures are 

7271 - 7i72 



7il7l' 



k\ = k 
The area is 



(meridian) 



72 



7i I 



(parallel) 



p r2ir rl rl rl 

|E|= dS= / ^/^{t)dtds = 2TT 7l (t)\j(t)\dt= / d/x 7 , (2.14) 
is Jo Jo Jo Jo 



and the volume enclosed by the surface is 



Vol(S) =vr [\ 2 (t)j 2 (t)dt. 
Jo 



It is then straightforward to check that Helfrich energy for axisymmetric surface and 
phase is given by expressions (1.10)-(1.12). 

2.3 Functions of bounded variation 

Let U denote an open subset of R; following [14] we say that a function / € ^(U) has 
bounded variation in U, and write / G BV(U), if 



supjy /y/cte : V € Cl(U), \ij>\ < l| < 00. 
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We say that a function / G L\ oc {U) has locally bounded variation in U, and write 
/ G BV\ oc (U), if for each open set V CC U 

supjy /^'dx : V € C£(J0,M < l| < oo. 

If / € BV\ oc (U) there exists a finite Radon measure // and a /x-measurable function 
77 : {/ ->■ R such that 



|r/| = 1 fi-a. e., and / fip'dx = — tpn 

Ju Ju 



dfx 



for all V G C£(V). We write \\Df\\ = n and [Df] = r]\\Df\\. If / G BV(t7)> then for each 

y c [/ 

P/IIOO = sup {J v f1>'dx : V e C^F), I VI < l} < 00. 

(Since D/ is scalar, we should have written \Df\ instead of but \Df\ and [Df] 

look too similar). For example, denoting by 5 X the Dirac distribution centered at x, if 
/ G BV(U) is the characteristic function of an interval (a, b) C U, then [Df] = S a — 5b , 
\\Df\\=S a + S b , 

[ ftfdx = - [ 1>d[Df\=Mb)-1>(a), [ ^d\\Df\\=i>(b)+i>{a) 
Ju Ju Ju 

and \\Df\\(U) =n°{a,b} = 2. 

In the sequel we will rely on the two fundamental results ([14, Section 5.2]) 

Theorem 2.2 (Lower Semicontinuity in BV). Let f k ,f G BV(U) be such that f k — > / 
inL\ QC {U). Then 

hmmi\\Df k \\(U)>\\Df\\(U). 

Theorem 2.3 (Compactness in BV). Let f k G BV(U) be such that snY> k {\\f k \\ L i{jj) + 
||-D/ fc || (U)} < 00. Then there exist a subsequence f k i and a function f G BV(U) such 
that f k > ^ f in L l {U). 

2.4 A notion of convergence for measure-function couples 

We recall that a sequence of Radon measures fj, n is said to converge weakly-* to fx G 
RM(R) if 

lim f <f>(t) dfi n (t) -»• f (f>(t) dfi(t) 

for every 4> G C C (M). We define the space of p-summable functions with respect to a 
positive Radon measure /j, as 

L p (n;M. 2 ) := j/ : R -»• M 2 ^/-measurable, such that ^ |/(x)| p ^(x) < +00 j . 
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Definition 2.2. [Convergence of measure-function couples] Following [2, Definition 5.4.3], 
given a sequence of measures fi n G RM(R) converging weakly-* to fx, we say that a 
sequence of (vector) functions f n G L 1 ( / u n ; R 2 ) converges weakly to a function / G 
L l (ji; R 2 ), and we write f n / in L^jR 2 ), provided 

lim / /" (*) • <K*) -> / /(*) • 0(t) d/i(t) (2.15) 

n ^°° Jr Jr 

for every <f> G C£°(IR;R 2 ). For p > 1, we say that a sequence of (vector) functions 
f n G F^";! 2 ) converges weakly to a function / G L P (^;M 2 ), and we write f n — 1 / in 
LP(/i™;IR 2 ), provided 

sup [ \f n (t)\ p dfi n (t) <+oo and / n ->> / in L 1 ^"; M 2 ). (2.16) 

nGN JR 

For p > 1, we say that a sequence of (vector) functions /" G L p (fj 1 n ; R 2 ) converges strongly 
to a function / G L^jM 2 ), and we write f n -»• / in #>(/x n ;R 2 ), if (2.16) holds and 

limsup||/ n || LP( ^. R2) < II/ILp^) ■ 

n— >oo 

Lemma 2.4 (Weak-strong convergence in L p (n;R d ) [20, Proposition 3.2]). Let p, q G 
(1, oo) such that 1/p + 1/q = 1. Suppose that fi n and fx are Radon measures on R and 
that f n G L P 0u n ;R d ) ; / G U>(n;R d ), g n G L% n ;M d ), g G L«(//;R d ) 6e swc/i that 

f n ^f weakly in L p (fi n ; R d ), g n -> 5 siron^ in L 9 (/i"; R d ). 

T/ien 

f n 9 n ^ fg weakly in L\fi n ;R d ). 

Theorem 2.5 (Lower semicontinuity [2, Theorem 5.4.4 - (ii)]). Let p > 1, let f n G 

L p (/i";IR 2 ) be a sequence converging weakly to a function f G L p (n;R 2 ) in the sense of 
Definition 2.2, then 

liminf / g(f n {t))d^(t)> [ g(f(t))dfx(t), 
n ^°° Jr Jr 

for every convex and lower semicontinuous function g : R — >■ (— oo, +oo]. 



3 Existence of a minimizer 

The proof of existence relies on the following fundamental estimate, which shows that, for 
properly chosen coefficients, the functional ^"(S, (p) bounds from above the L 2 -norm of 
the second fundamental form of S, independently of the phase <p. Whereas the estimate 
follows from the phase-independent case (see [7, Lemma 2.1]), for sake of completeness 
we include the details. 

Lemma 3.1 (Fundamental estimate). If k 1 h ,k 1 g , i = A,B satisfy 

4/>0 and ^G(-2,0) fori = A,B, (3.1) 
K H 
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then there exists C > such that 



J^(E, f)>C^J {kj + k 2 2 )dS - |E|^ 



for all admissible couples (E, <p). 
Proof. Let A, Ho G R, note that 

+ fc 2 ) 2 + Afcife = i(fc? + fc£) + (1 + A)fcifc 2 > l ~\\ + X \ k{ + fc 2 ), 

and the last term is positive if and only if A € (—2, 0). For all e > it holds 

H 2 {H _ Ho + Ho) 2 1+ l + £ 

~2~ ~ 2 <— (^-^o) +^^0, 

and thus 

i±V - «,)» + ^ + A(l + s)K > 1 - |1+ 2 A(1+e)l ( t ? + «). 

For all k// > 0, choosing A = kq/hh S (—2,0) and e > such that (1 + e)kg/kh G 
(-2,0), we get 

^ (tf - tf ) 2 + kqK + Cl tf 2 > c 2 (A; 2 + fc 2 ), 

where ci = Kif/2e and c 2 = Kg ~ l ^ 1 + _g (1+e)l > 0. Denoting by c\,c\ (i = A, B) the 
constants obtained for kh = k % u , kg = k 1 g , etc., we integrate on E to obtain 



^(E, = jf f^^-(H - H (<p)) 2 + K G (p) + ^(r) 

> min {c^, cf } ^(fc 2 + k 2 2 )dS - |E| max {cf(H^) 2 , cf (i7 B ) 2 } 



(3.2) 



□ 



3.1 Compactness 

Let the total area, the ^4-phase area and the volume constraints A, ILa, V be given. Let 

S n = |(E^,^),...,(S" (n) ,^ m(n)) )| , 
where m(n) is the cardinality of the rt-th family, be a minimizing sequence for JP, i.e., 
S n G S(A, Ha, V) and liminf F{S n ) = inf {F{S) : 5 € <S(.A, U A , V)} . 

n— >co 
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Compactness for curves. Since for any ^4, IIa, V satisfying (1.13) it is possible to 
construct a spheroid with area A, volume V, and divide its surface in two domains such 
that one has area IL4, the infimum above is finite, and there exists A > such that 
T(S n ) < A for all n G N. In [7] we studied the compactness properties of a sequence of 
surfaces for which the bound 

J (kf + kf) dS < C (3.3) 

holds uniformly. Owing to Lemma 3.1, we can apply the results in [7, Lemma 2.4, Lemma 
2.5, Lemma 3.6, Proposition 3.7] to the sequence of surfaces {S^}. We summarize these 
results in the following proposition. Regarding notation, we denote the first and the 
second component of a curve jj by (77)1 and (7^)2 • 

Proposition 3.2. Let S n be a sequence of finite systems of admissible surfaces and 
phases satisfying the bound (3.3) for some fixed constant C > 0. Denote the generating 
curves by (7™, . . . , 7^/ n \)- Then, there exists a subsequence S nk such that 

(i) the cardinality m(n-fc) of the system is uniformly bounded. Therefore, it is not 
restrictive to assume that m(nk) = cj, for some constant u > 0. 

(ii) There exists a system of curves 71, . . . , 7j such that for all j = 1, . . . , J (J < u) 

7 j eW 1 > oo ((0,l);M 2 h ), ( 7j )i>0 a.e. in (0,1), ^ e L 2 (^;R 2 + ), 

and, up to a permutation of the indices, 

7™ fe 7j strongly in H l ((0, 1);R+), (3.4) 

7™ fe Jj weakly in L 2 (p,^ k ;^\). (3.5) 

(Hi) For j = J + 1, . . . ,lo, the curves 7™ fe converge to a point lying on the z-axis, i.e. 

(7™ fc )i -> 0, (77* ) 2 -> ^ , sfron^/y in ^((0, 1)) 
/or some Zj £l. 

For j = 1, . . . , J, it holds = 0} < +00 and 7j G W 2,2 ({(7j)i > 0};K 2 ) 

(vj For j = 1, . . . , J, each curve jj meets the z-axis orthogonally, so, in particular, the 
generated surface £j is the union of a finite number of C 1 -regular surfaces. 

(vi) The area and volume constraints pass to the limit, i.e., 

m(n) J m(n) J 

E W\ = Ei^i> E voi(^) = Evoi(^). 

j=i j=i j=i j=i 

Since singularities for jj can occur only on the z-axis, i.e. where (7^)1 = 0, points (iv) 
and (v) imply that the limit system can be reparametrized as a finite family of admissible 
curves [7, Corollary 2.9]. 
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Compactness for phases. We turn now to the question of compactness for the phases. 
Since the system S n has fixed cardinality, it is enough to study the behavior of a single 
couple. Let (Y, n ,ip n ) be a sequence of admissible surface-phase couples, where S n is 
generated by 7™, and let 7 be a curve such that 7™ — > 7 as in Proposition 3.2 (ii), (iv). 
First of all, since ip n takes values in {0, 1}, we can find a subsequence (not relabeled) 
and a function <p G L°°(0, 1) such that 

<p n ^ ip weakly-* in L°°(0, 1), < ip(t) < 1 for a.e. t G (0, 1). (3.6) 

Since convergence (3.4) implies that 7™ — > 71 uniformly in C°([0, 1]), for every compact 
K C {71 > 0} there exists e > such that 7™ (t) > e for all t G K. More precisely, since 
by Proposition 3.2-(iv) the set {71 = 0} is finite, for every 5 > there exists a compact 
set K = K(5) C {71 > 0} and a positive number e = e(S) > such that 

£^[0, 1]\K) < - and 7l n (t) > e VtGif. (3.7) 
Owing to Lemma 3.1 

A > J^(£ n , p n )>C^J (k[ 1 ) 2 + (A£) 2 dS - |£ n f) + aT^P 1 ). 

Recalling (1.6) we obtain 

A ^ 4 >^ 1 (r")= Yl ^(0> E 7rW>^°(J(^)nif). 

teJ(ip n ) teJ{ip n )c\K 
Let U C K be an open set such that £ 1 (^\C7) < 5/2, we have that ip n £ BV(U) and 

IKW) = \\<P n \\LHu) + \\D<p n \\(U) 

= / p n (t)dt + n°{J(^ n )nU) < 1 + 



27ra£ 

By classical compactness and lower semicontinuity for BV functions (see Theorems 2.2 
and 2.3), there exists a subsequence ip nk and a function (p G 5V(?7) such that 

cp nk ^p strongly in ^(U), (3.8) 
liminf ||D(p" fc ||(?7) > ||D^||(17). (3.9) 

By (3.6), we also have that <p = tp a.e. in U and thus is the strong limit for the whole 
sequence ip n . By (3.7) and (3.8), we found that for every 5 > there exists an open set 
U C (0, 1) and a function (p such that 

lim f \<p n -ip\dt= lim ( f \<p n -^\dt+ [ \<p n -v\dt) 
n ^°° Jo n ^°° \Ju J[o,i]\u J 

^ n -Cp\dt + 25 



= lim 1 


u 


n^oo 




< lim I 




n^oo 





< 25. 
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Since 5 is arbitrary, we obtain strong convergence in L 1 (0, 1). Recalling that strong 
convergence, up to subsequences, implies convergence almost everywhere, we conclude 
that convergence (3.6) is improved to 

^ n -> <p strongly in L\0, 1), 99 G W loc ({ 7l > 0}; {0, 1}), 

or, actually, noting that \ip n \ < 1, to 

<p n -> v? strongly in L p (0, 1), y> G BV loc {{ 7l > 0}; {0, 1}), 

for all p G [l,oo). Note that ^-convergence for the phases, combined with convergence 
(3.4) for the curves, implies that the phase area constraint passes to the limit: if S n = 

i=i J ° i=i J ° 

and since ( 7 i)i = for i = J + 1, . . . , u, 

J 



V2vr / <p i (t)(i i ) 1 (t)\ji\dt = n A . 
i=i Jo 



We collect these computations and the results of Proposition 3.2 in the following 

Proposition 3.3 (Compactness). Let the area, A-phase area and volume constraints 
A,Ha,V be given. Let 



^ = {(S™^),...,(S« (n) ,^ m(n)) )} 



be a sequence of systems of admissible surfaces-phases such that S n 6 S(A, 11a, V) and 
^(S 71 ) < A, for some constant A > 0. Then there exists a subsequence S nk and an 
admissible system 

S = {(S l5 ^), . . . , (Ej, ipj)} G 5(A n A , V) 
such that for all j = 1, . . . , J 

VeW 1 ' 00 ^);®*.), hj)i>0 a.e. in (0,1), 7j G L 2 (/x 7 ,;M^), 

^G W loc ({( 7j )i>0};{0,l}) 
and, up to a permutation of the indices, 

7™ fe -»• 7,- strongly in ^((0, 1);R+), (3.10) 

7?* ~^ 7j weakly in L 2 (n n fc ;R+), (3-H) 



'3 
and 



<p™ k -»• strongly in L p (0, 1), Vp G [1, 00). (3.12) 



liminf ||Z?VJ nfc || (i7) > H^ll(^) Vf/CC{7i>0}. (3.13) 
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3.2 Lower semicontinuity 

In this Section we prove the following 

Proposition 3.4 (Lower semicontinuity). Let (T, n ,p n ) be a sequence of admissible cou- 
ples of surfaces-phases, converging to the couple (S, ip) in the sense of (3.10)-(3.12) ; 
then 

liminf^"(£ n ,^ n ) > Jf(S,y>). 

n— >-oo 

Proof. We first address the lower semicontinuity of the term aH 1 ^). We need to prove 
that if the couple (-y n ,ip n ) converges to (7, <p) as in (3.10)-(3.12), then 

liminf aft 1 (r") = liminfc727r [ 7? d\\D<p n \\ > a2n [ 7l d\\D<p\\ = a^CT). (3.14) 

rwoc rc->oc J Q J Q 

This result would be straightforward, if we had ip n <p in BV(0, 1). Since p n is only 
BV\ oc {{li > {0' 1})' we need one remark: 

Lemma 3.5. Let £1 C R 6e open and bounded. Let 7,7™ G Co(^) 6e nonnegative, and 
ip, ip n £ i?V5 oc (f2) oe swc/i 

7™ ->■ 7 uniformly in C°(0), (3.15) 

liminf ||^ n ||([/) > H-D^ll(^) forallUCCQ. (3.16) 

n— >oo 

7/ £/tere exists C > swc/i i/iai 

sup / 7 n d||ZV*ll < C, (3.17) 
« isi 

then ^WDifW is a finite Radon measure on and 

liminf [ j n d\\D<p n \\> f jd\\D<p\\. 
Jn Jn 

Proof. Let i\) G C c (0), with < V < 1, 

/ 7 n d||ZV 1 ||> f Y l i>d\\Dip n \\= [ ( 7 n - 7 )^d\\Dip n \\+ [ jipd\\D<p n \\. 
Jn Jn Jn Jn 

Since p n £ BV\ oc (Q), there exists M = M(ip) such that 

/ ipd\\D<p n \\ < M. 
Jn 

By (3.15), for all e > there exists n such that sup^ (7™ — 7I < e for n > fi. We have 

liminf [ -f n d\\Dip n \\ > liminf f- / | 7 n - 7^ d||£V n || + / T^dp^" 
n -*°° n ^°° V in Jn 

(3.16) /• 

> -Me+ / 74>d\\Dip\\. 
Jn 
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Since e was arbitrary, we obtain 

liminf [ -j n d\\D^ n \\> f "/il>d\\D(p\\, G C c (0), < < 1. 

n ^°° Jn Jn 

Let now £ C c (f2) be a sequence such that ^ k (x) = 1 if x) > < ^(x) < 1. 

Since 7 6 Cb(fi), there exists a continuous and monotone increasing function w : — >■ M 
such that w(0) =0 and 

sup\j(x)ip k (x) -j{x)\ < sup {7(0;) : d(x,dQ) < - 1 < w(l/fc), 

i.e., ^ip k —7- 7 uniformly in Cq{VL). Since the mapping 1— > J*^ 7^ d||D^|| is continuous 
with respect to uniform convergence, we can extend 7||-D<^|| to a finite Radon measure 
on Q and conclude that 

f -fd\\Dip\\= lim f ji> k d\\D<p\\ < liminf f j n d\\D(p n \\. 

Jn k^ooj n n->oo J n 

□ 

By (3.10)-(3.13) we can apply Lemma 3.5 to a sequence of admissible couples and 
obtain the lower semicontinuity (3.14). 

Regarding the remaining part of the functional, define 
(£, if) : = ±J K H {y){H - H {cp)) 2 dS + J KG (<p)KdS 

= ir k h (<p) {h + k 2 - Ho(ip)) 2 -/tljldt + 2ir kg(<p)&i&27i|7| dt. 
Jo Jo 

We need to show that if the couple (7™, ip n ) converges to (7, tp) as in (3.10)-(3.12), then 

liminf ^(S ra ,v? n ) > J4?(Z,<p). (3.18) 

Define the Radon measures 

li-yn ■= 2iry? \j n \ £^[0,1] , ^7 := 2tt7i | 7 | £^[0,1] , 

A" := irK H { V n )^\T\C\ m , r := 2™ G ( ( ^hr|7 n |£ 1 L [0 , 1] ■ 

By (3.12) and by the linearity of kh , kg , -f^o , 

K H (ip n ) -»• K H (<p), kg(<P h ) -> «*?(¥>), #0(9?™) -> #0(9?) 
strongly in L p (0, 1) for every p G [1, +00), while by (3.10), 

7i n |7 n | -> 7i|7l uniformly in C°([0, 1]), 

so, in particular, A n — ^ A := |/ejj((p)7i|7|jC 1 L[o j i] and tp n ^ ip := kg-(^)7 1 |7|£ 1 l[o i i] as 
measures. Recalling also (3.1), it is straightforward to check that weak (resp. strong) 
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convergence in L p (/i 7 n) is equivalent to weak (resp. strong) convergence in L p (X n ), or in 
L p (ip n ), in the sense of Definition 2.2. By convergence (3.11) and [7, Lemma 3.4], 

H n — H weakly in L 2 {X n ), K n -» K weakly in L 1 ^™). 

By Theorem 2.5 and Lemma 2.4 we conclude that 

liminf f — (cp n )(H n - H (<p n )) 2 d^n + [ K G (cp n )K n d^n 

n— >oo J 2 J 

= liminf f (H n - H (<p n )) 2 d\ n + f K n d^ n 

n^oo J J 

> J(H - H (p)) 2 d\ + J Kdij> 

= j -^{^>){H - H ) 2 d^ + J K G ((p)Kd^. 
This proves (3.18), which together with (3.14) concludes the proof of Proposition 3.4. □ 
3.3 Proof of Theorem 1.1 

Let the total area, A-phase area and volume constraints A, n^,V be given, such that 
the isoperimetric inequality and the bound on the phase (1.13) are satisfied. Let the 
parameters n l H , k g satisfy (1.9) and H'q £ R, for i = A, B. Let the set S(A, ILa, V) and 
the functional F be given as in the statement of Theorem 1.1. 

Let S n = ((£?, ip n ), ... , (S™ (n) , ip n m[n) )) £ S(A, U A , V) be a sequence of admissible 
couples of surfaces and phases such that 

lim F(S n ) = inf ?(S). (3.19) 

n^oo S&S(A,U A ,V) 

Since S n satisfies J 7 (S n ) < A, for a suitable A > 0, by Lemma 3.1 there is a constant 
C > such that 

We can therefore apply Proposition 3.3, and find admissible couples (71, ipi), . . . , (77, ipj) 
and a subsequence (not relabeled) of constant cardinality m{n) = m, such that for 
< j < J, as n — > 00 

7j -> 7j , -> ¥>j 

in the sense of convergence (3.10), (3.12), while 7j +1 , • ■ . ,7^ shrink to points, thus not 
contributing to the total area of the system. The system S of surfaces-phases couples gen- 
erated by ((71, ipi), ■ ■ ■ , (7j, <£j)) satisfies the total area, phase area and enclosed volume 
constraints. By the lower semicontinuity Proposition 3.4 

liminf F(S n ) > F(S), 

n->oo 

so that, by (3.19), T(S) = inf J 7 . The proof of Theorem 1.1 is thus complete. 
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